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Abstract. The reionization of the Universe is expected to 
leave a signal in the form of a sharp step in the spectrum 
of the sky. If reionization occurs at 5 ^ z- lon £ 20, a feature 
should appear in the radio sky at 70 v ^ 240 MHz due 
to redshifted Hi 21-cm line emission, accompanied by an- 
other feature in the optical/near-IR at 0.7 A 2.6 /im 
due to hydrogen recombination radiation. The expected 
amplitude is well above fundamental detection limits, and 
the sharpness of the feature may make it distinguishable 
from variations due to terrestrial, galactic and extragalac- 
tic foregrounds. 

Because this is essentially a continuum measurement of 
a signal which occurs over the whole sky, relatively small 
telescopes may suffice for detection in the radio. In the 
optical/near-IR, a space telescope is needed with the low- 
est possible background conditions, since the experiment 
will be severely background-limited. 
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early Universe - diffuse radiation - observations 



1. Introduction 

The epoch of reionization marked the end of the "dark 
ages" during which the ever-fading primordial background 
radiation cooled below 3000 K and shifted first into the 
infrared and then into the radio. Darkness persisted un- 
til early structures collapsed and cooled, forming the first 
stars and quasars that lit the universe up again. During 
this epoch the volume filling factor of ionized hydrogen 
(Hn) increased rapidly. It is a generic feature of theo- 
retical models and three-dimensional numerical simula- 
tions that reionization occurred at 5 ^ z lon ^ 20 in most 
cold dark matter (CDM) cosmogonies (Gnedin & Ostriker 
1997; Haehnelt & Steinmetz 1998; Cen 1998; Haiman & 



Loeb 1999) Q We know reionization took place before 
z ss 5, as evidenced by the lack of hydrogen continuum ab- 
sorption in the spectra of high-redshift quasars (Schneider 
et al. 1991) and galaxies (Franx et al. 1997). It is unlikely 
that reionization occurred at z ^ 50, for in that case the 
level of degree scale anisotropy in the CMB would be lower 
than observed on ten degree scales (e.g. Knox et al. 1998). 
Tilted and neutrino dominated CDM models give z lon ~5- 
8, and open and A-dominated models give z- mn ~7-30 (Cen 
1998). 

Moreover, since the characteristic distance between the 
sources that ionized a rapidly recombining intergalactic 
medium (IGM) was much smaller than the Hubble radius, 
the transition from H I to H n was quite abrupt; the over- 
lapping of the H II regions surrounding the first star clus- 
ters and miniquasars was completed on a timescale much 
shorter than the Hubble time. The reionization epoch oc- 
curred almost as a "phase transition" of the Universe. 
From the observational point of view, this poorly con- 
strained epoch is among the most important unknown 
quantities in cosmology. 

Various techniques for probing the history of the tran- 
sition from a neutral IGM to one that is almost fully ion- 
ized have been proposed in the literature. One way is to 
look for the 21-cm hyperfine line of neutral hydrogen, red- 
shifted to frequencies in the range 70-240 MHz (Madau et 
al. 1997; Gnedin & Ostriker 1997; see also Swarup 1984; 
Swarup & Subrahmanyan 1987; Scott & Rees 1990). Prior 
to the reionization epoch, the neutral gas that had not yet 
been engulfed by an H n region may be seen as 21-cm line 
emission if a mechanism existed that decoupled the spin 
temperature from the CMB temperature. Physical mech- 
anisms that would produce a 21-cm signature are Lya 
coupling of the spin temperature to the kinetic tempera- 
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Note, however, that while simulations are able to track the 
formation and merging of dark matter halos and the subse- 
quent baryonic infall, they are much less able to predict the 
efficiency and rate of radiation emission from gravitationally 
collapsed objects. 
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ture of the gas, preheating by soft X-rays from collapsing 
dark matter halos, and preheating by ambient Lya pho- 
tons (Madau ct al. 1997). 

So far the emphasis has been on detecting fluctuations 
due to individual large Hi concentrations. This is very 
difficult in the presence of a strong fluctuating foreground 
emission; the H I lines from individual concentrations must 
be resolved, and this requires large telescopes with high 
sensitivity and spectral and angular resolution. Alterna- 
tively, if reionization occurred abruptly and the fraction 
of neutral hydrogen underwent a drop of a factor of 10 3 in 
about a tenth of a Hubble time as predicted (e.g. Gnedin 
& Ostriker 1997; Baltz et al. 1998), it is conceivable that a 
global signature may be detectable with telescopes of mod- 
est size in the extragalactic background spectrum at meter 
wavelengths. This is in principle a much simpler measure- 
ment. Some existing ground-based telescopes may be suit- 
able, or, if terrestrial interference is too severe, a dedicated 
space project might be considered. A complementary tech- 
nique is to look for a signal in the optical/near-IR back- 
ground due to hydrogen recombination radiation from the 
reionization epoch (Baltz et al. 1998). 

The prospects for detecting a global signature of the 
reionization epoch are considered below. In the following 
we shall denote the present-day Hubble constant as H = 
100 h km s- 1 Mpc" 1 . 

2. Reionization signature 

Numerical N-body/hydrodynamic simulations of struc- 
ture formation in a strongly clumped IGM have started 
to provide a picture for the origin of intervening absorp- 
tion systems, one of an interconnected network of sheets 
and filaments, with virialized systems located at their 
points of intersection. The gas clumping factor rose above 
unity when the collapsed fraction of baryons became non- 
negligible, i.e. at z ^ 20, and grew to a few tens at 
z « 8 (Gnedin & Ostriker 1997). This made the volume- 
averaged gas recombination time, 

t rcc = [1.17n H a B C]" 1 = 0.08 Gyr (j^j 

shorter than that for a uniform IGM and shorter than the 
the Hubble time at that epoch. In the above expression, 
nH is the mean hydrogen density of the expanding IGM, 
n H (0) = 1.7 x 10~ 7 (Ct b h 2 /0.02) cm~ 3 , a B is the recom- 
bination coefficient to the excited states of hydrogen, and 
C = (^hiiV^hii i s t ne ionized hydrogen clumping factor. 

When an isolated point source of ionizing radiation 
turns on, the ionized volume initially grows in size at a rate 
fixed by the emission of UV photons, and an ionization 
front separating the Hn and Hi regions propagates into 



the neutral gas. The evolution of an expanding H II region 
is governed by the equation 

(2) 

dt flH Tree 

(Shapiro & Giroux 1987), where Vj is the proper volume 
of the ionized zone, iVj on is the number of ionizing photons 
emitted by the central source per unit time, and H is the 
Hubble constant. Across the ionization front the degree of 
ionization changes sharply on a distance of the order of the 
mean free path of an ionizing photon. When t rcc <C t, the 
growth of the H n region is slowed down by recombinations 
in the highly inhomogeneous IGM, and its evolution can 
be decoupled from the expansion of the universe. As in 
the static case, the ionized bubble will fill its time- varying 
Stromgren sphere after a few recombination timescales, 

Vi = N ion t lec e _ t/i 
n H 

While the volume that is actually ionized depends on the 
luminosity of the central source, the time it takes to pro- 
duce an ionization-bounded region is only a function of 

^rcc • 

At early times, the characteristic distance between the 
sources that ionized a rapidly recombining IGM was much 
smaller than the Hubble radius. The phase transition from 
Hi to Hll was therefore quite abrupt. When the Hn re- 
gions surrounding the first star clusters and miniquasars 
overlapped, the fraction of neutral hydrogen dropped by 
a few orders of magnitude in about a tenth of the Hubble 
time at that epoch, the continuum optical depth decreased 
suddenly and the background ionizing intensity underwent 
a drastic increase in a very small redshift interval (Gnedin 
k Ostriker 1997; Baltz et al. 1998). The redshift of this 
event is expected to be essentially the same in all direc- 
tions. It is this abrupt change in Hi absorption that flags 
the epoch of reionization and may be observable at radio 
and optical/near-IR frequencies, as shown below. 

3. Detectability of an H I edge 

The intergalactic medium prior to the epoch of full reion- 
ization should be detectable in 21-cm line radiation. In 
the absence of a decoupling mechanism, the spin temper- 
ature of neutral hydrogen would go to equilibrium with 
the CMB, and no emission or absorption relative to the 
CMB would be detected. The collapse and cooling of the 
first non-linear structures had a twofold effect: (1) it pre- 
heated the IGM to temperatures ~ 150 K, above that of 
the CMB; and (2) through scattering by Lya photons it 
provided a mechanism to couple the spin temperature to 
the kinetic temperature of the gas (Madau et al. 1997). A 
patchwork - due to large-scale structure and non-uniform 
heating and coupling - of 21-cm line emission would result, 
a signal which would have disappeared after reionization. 
A signature ("step") in the continuum spectrum of the 
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radio sky, at a frequency of ~ 70 — 240 MHz for z- lon in 
the range 5-20, will therefore flag the reionization epoch. 

To illustrate the basic principle of the observations we 
propose, consider a neutral IGM with spin temperature 
Ts ^> Tqmb- In an Einstein-de Sitter universe (Oo = 1> 
= 0), its intergalactic optical depth at 21(1 + z) cm 
along the line of sight, 

^^-(Z^) («£)(! + ,)•/>, ,4) 

will typically be much less than unity. The experiment 
envisaged consists in its simplest form of two measure- 
ments separated in frequency such that the one at shorter 
wavelengths detects no line feature because all hydrogen 
has been reionized. The differential antenna temperature 
across the "step" , as observed at Earth, will be 

AT = (1 + zY^Ts - T CMB )(1 - e- T ) 

where Zi on is the redshift of the transition epoch. Note 
that in a universe with Qq -C 1, AT increases more nearly 
linearly with (1 + z ion ), and the numerical coefficient in 
Eq. (5) may be larger by up to a factor 3[(1 + Zi on )/9] 1 / 2 . 
Depending on the cosmology, AT will increase toward low 
frequencies at a rate between v~ x l 2 and v~ x . Hence, the 
expected cosmological signal will be an edge of more than 
~ 0.02 K (for h = 0.5) at 70-240 MHz superimposed on 
the 2.73 K cosmic background radiation. Fig. 1 shows the 
total cosmological signal, comprised of the cosmic back- 
ground radiation with the reionization step superimposed. 

3.1. Fundamental sensitivity limits 

The fundamental limitation on the sensitivity achievable is 
set by the intensity of the galactic and extragalactic fore- 
grounds. Modern receivers are such that the system tem- 
perature should be largely dominated by these extrater- 
restrial signals. The limiting sensitivity is thus 

AT min ~ -7=, (6) 

where T sys is the total system temperature, 8v is the band- 
width, and t is the integration time. For T sys = 150K (the 
coldest regions of the sky at 150 MHz) with a bandwidth 
of 5 MHz and integration time 24 hours, AT mm <~ 0.0002 
K. The reionization step predicted by Eq. (5) would be 
~ lOOer, independent of telescope size. In a real measure- 
ment the system temperature would of course be higher 
than this fundamental limit of 150K, but clearly sensi- 
tivity is not an issue, and the challenge concerns signal 
contamination and calibration, as discussed below. 
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Fig. 1. Expected brightness temperature of the cosmic back- 
ground in the vicinity of the Hi reionization edge as a func- 
tion of observing frequency. Three cases are shown for the Hi 
step, (a): the case from Eq. (5) with z ion = 7, Az — 1 and 
h — 0.7, (b) and (c): revised results provided by N. Gnedin 
(private communication) from the simulations by Gnedin & 
Ostriker (1997) and Baltz et al. (1998) case A respectively, 
with n M = 0.35, ft a = 0.65, and h = 0.7. 



3.2. Contamination by galactic and extragalactic 
foregrounds 

The difficulty posed by the galactic and extragalactic fore- 
grounds is that they can be complex, both in frequency 
and position. The question is whether one can in princi- 
ple extract a 0.02K step from this much stronger varying 
continuum. 

All-sky maps at 150 and 408 MHz are presented in Lan- 
decker & Wielebinski (1970) and Haslam et al. (1982). The 
regions of lowest brightness temperatures occur in large, 
relatively smooth regions away from the galactic plane 
and galactic loops. This emission is comprised of four 
components: galactic synchrotron emission (~ 70% at 150 
MHz), galactic thermal (free-free) emission (~ 1%), the 
integrated emission from extragalactic sources (~ 27%), 
and the 2.73 K cosmic background itself. 

The galactic nonthermal synchrotron emission is dom- 
inant at these low frequencies, and has been studied since 
the pioneering days of radio astronomy. Its spectrum is 
close to a featureless power law, although there are grad- 
ual variations in the spectral index with position on the 
sky and with frequency. The spectral index is smallest in 
the coldest regions of the sky, away from the galactic plane 
and galactic loops. Variations appear to be mainly due to 
the latter. At 100 MHz the spectral index (3 of this com- 
ponent is about -2.55 in cold sky regions, and it steepens 
to -2.8 at 1 GHz (Bridle 1967; Sironi 1974; Webster 1974; 
Cane 1979; Lawson et al. 1987; Reich & Reich 1988; Ban- 
day & Wolfendale 1991; Platania et al. 1998). In some 
regions such as the north galactic pole this steepening of 
the spectrum appears to be considerably less pronounced 
(Bridle 1967). 
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There may also be a dispersion in the spectral index 
at each position on the sky, due to distinct components 
along the line of sight. Some indication of this dispersion 
is given by point-to-point variations of the spectral index 
over the sky (Landecker & Wielebinski 1970; Milogradov- 
Turin 1974; Lawson ct al. 1987; Reich & Reich 1988; Ban- 
day & Wolfendale 1990, 1991). It is not large, and appears 
to be smallest at the lower frequencies (Lawson et al. 1987; 
Banday & Wolfendale 1990). A value of <t(j3) ~ 0.1 is prob- 
ably appropriate at 100-200 MHz. Further information on 
small-scale uniformity can be provided by simultaneous 
interferometric data on many short baselines using tele- 
scope arrays such as the WSRT and GMRT. The recently 
completed Westerbork 325 MHz survey of the northern 
hemisphere (WENSS, de Bruyn et al. 1998) already gives 
some constraints. At high galactic latitudes, where the 
minimum brightness temperature of the diffuse galactic 
emission is about 20 K at 325 MHz, the variations in the 
brightness temperature of the galactic foreground emission 
are typically less than about 0.4 K, or about 2%, on scales 
ranging from 5-30 arcmin. If spatial intensity variations 
are accompanied by spectral variations (as they would be, 
if caused by low surface brightness H II regions or super- 
nova remnants), such measurements can be used to set 
limits on spectral variations. 

The galactic thermal emission at high galactic lati- 
tudes is due to very diffuse ionized gas, with total emission 
measure about 5 pc cur 6 for T e = 8000K (Reynolds 1990). 
More recent studies of the properties of this gas have been 
made by Kogut et al. (1996) and de Oliviera-Costa et al. 
(1997), who show that it has structure correlated with 
high-latitude dust clouds. This gas is optically thin at fre- 
quencies above a few MHz, so its spectrum in the region 
100-200 MHz is well-determined, with a spectral index of 
-2.1. 

Spectral lines present a related potential contamina- 
tion. Galactic radio recombination lines occur every 1-2 
MHz over the frequency range of interest. In addition to 
spontaneous emission from ionized hydrogen and other el- 
ements in the clouds, there is also stimulated emission (or 
absorption) against the nonthermal background. Peak line 
intensities can reach IK or so, although it turns out that 
the galactic ridge recombination lines happen to make the 
transition from emission to absorption at 100-200 MHz, 
and so are at a minimum in this range (Payne et al. 1989). 
Even a line intensity of IK becomes diluted to 0.002K in 
a 5 MHz band. In any case, spectral lines can be iden- 
tified and removed with observations of higher spectral 
resolution. Observations such as these may provide the 
added bonus of discovering or setting limits on presently 
unknown spectral lines. 

The isotropic emission from extragalactic sources has 
been estimated both directly and from integrated source 
counts. Bridle (1967) found that this component amounts 
to about 48K at 150 MHz, with a spectral index of -2.75. 
Cane (1979) obtained a value corresponding to about 32K 



at 150 MHz, based on observations at 10 MHz. This is 
consistent with results from integrating the source count 
data (Willis et al. 1977; Simon 1977; Lawson et al. 1987). 
Simon (1977) estimated that the total extragalactic back- 
ground should turn over between 1 and 2 MHz due to syn- 
chrotron self-absorption in the individual source compo- 
nents, but individual sources can have complex frequency 
structure at much higher frequencies. In the extreme, spec- 
tral indices can range from +2 (sources with sharp syn- 
chrotron self- absorption turnovers) to -3 (pulsars). Over- 
all, the resulting features in the extragalactic background 
should be small and slowly varying functions of frequency 
due to source inhomogeneities and the average over many 
sources. 

This has been simulated using the known spectral 
properties of 3C sources as studied by Kellermann et al. 
(1969). They classified the spectra of almost all 3C sources 
in the following categories: S (straight - single power law), 
C- (convex - steepening at higher frequencies), C+ (con- 
cave - flattening at higher frequencies) , and Cpx (complex 
- containing multiple components). Of the 299 extragalac- 
tic sources studied, 42% are S, 52% are C-, 3% are C+, 
and 3% are Cpx. The actual spectral shapes of 37 3C 
sources representing all these classes were measured from 
spectra published in Kellermann (1966, 1974), Kellermann 
et al. (1969), Kellermann & Pauliny-Toth (1969), and 
Kellermann & Owen (1988). These were then weighted in 
accordance with the frequency of occurance of the various 
classes, and co-added to give a representative composite 
spectrum of the extragalactic foreground. 

The resulting extragalactic spectrum in the range 50- 
300 MHz is close to a power law of spectral index -2.65 
(the index is steeper at higher frequencies). The deviation 
from this power law is a slowly and continuously varying 
function of frequency, with an index ranging from -2.59 
to -2.66 and no sharp features. The smoothness of this 
spectrum based on just 37 sources (see Fig. 3) is reas- 
suring, as vastly more sources will contribute to the real 
extragalactic signal, which will be correspondingly more 
featureless. The overall spectrum of the real extragalac- 
tic signal will differ from that computed here, as only the 
most luminous sources are represented in this simulation, 
but the important conclusion is that the spectrum should 
be smooth and featureless. 

Nevertheless it will still be important to identify 
and characterize the spectral properties of the dominant 
sources within the main beam. Interferometer observa- 
tions can be made for this purpose. A single source with a 
flux density of 35 mJy at 150 MHz would produce an an- 
tenna temperature increase of 0.01 K for a 45m telescope 
with 50% main beam efficiency. All sources should there- 
fore be identified down to at least this level. The GMRT 
would be the instrument of choice to conduct these high 
resolution observations. With receivers for 50, 150, 235 
and 325 MHz and a resolution of better than 1' (at 50 
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Fig. 2. Spectral index vs. observing frequency, as deduced over 
10 MHz intervals from simulated spectra which include the 
galactic nonthermal emission with median spectral index /3 = 
-2.55 and a Gaussian dispersion around this value with a = 
0.1, and the galactic thermal emission with EM = 5 pc cm 
and T e = 6000K, in addition to the reionization signal (the 
constant 2.73K cosmic background has been removed). The 
emission from extragalactic sources is omitted in this case. The 
total brightness temperature at 150 MHz is 150K. The points 
represent the case without a reionization step, and the error 
bars correspond to ±5cr for a 24-hour observation. The lines 
represent the three cases shown in Fig. 1. 

MHz), spectra of all significant sources can be determined 
in the relevant frequency range. 

Thus, it may be possible to quantify all the foreground 
contaminants sufficiently to model them and possibly re- 
move some of them. The fact that the reionization signal 
should be the same over the whole sky is a great advantage 
- searches in different positions on the sky with different 
contamination should give the same result. A range of 
different measurements can help greatly in disentangling 
these foreground signals. Accurate broadband spectra over 
a wide frequency range can help to determine the spectral 
behavior of the galactic nonthermal component, and mea- 
surements in adjacent positions can constrain the disper- 
sion of the spectral index. High frequency measurements 
can be used to accurately determine the galactic thermal 
component. Interferometer observations can be used to 
identify regions devoid of strong extragalactic sources, and 
to measure the properties of the most prominent sources 
in the regions chosen. 

The analysis can be done in a number of ways, de- 
pending on how effectively the different components can 
be modeled or removed. If the overall spectrum could be 
well determined at frequencies above that of the reioniza- 
tion step (e.g. v > 237 MHz, corresponding to z < 5), a 
simple extrapolation to lower frequencies could be made 
to find the step. It is unlikely, however, that such an ex- 
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Fig. 3. Spectral index vs. observing frequency as in Fig. 2, 
except that in this case the emission from extragalactic sources 
is included, represented by the simulation described in the text 
based on the actual spectra of 37 3C sources. Again, the total 
brightness temperature at 150 MHz is 150K. 

trapolation could be made with sufficient reliability. Alter- 
natively, after those components that are well-determined 
have been removed (at least the 2.73 K cosmic background 
and the galactic thermal emission), a best fit (power law 
or low-order polynomial) could be made and subtracted 
or divided from the actual spectrum to reveal the step, al- 
though the fit may well introduce artifacts that can mask 
or confuse the signal. A natural approach is to make the fit 
at the two ends of the spectrum and interpolate, but this 
requires a well-behaved spectrum and assumes that the 
signal is located near the middle of the spectrum. Another 
possibility is to measure the spectral index point-by-point 
along the spectrum, and look for an abrupt jump due to 
the step. 

Figure 2 shows the spectral index measured every 10 
MHz from simulated spectra containing just the galactic 
thermal and nonthermal emission in addition to the reion- 
ization signal. In Fig. 3 the extragalactic sources have been 
included. These two examples give an idea of the range of 
spectral curvature that can be expected. The overall cur- 
vature is substantial, but the reionization signal can still 
be seen in all three cases because of its relative sharpness. 

Another technique that does not involve any fitting 
procedure (but does require high-precision measurements) 
is trend analysis - analyzing the spectrum for deviations 
from a smooth trend. Fig. 4 shows the "trend ratio", de- 
fined here as the ratio of the observed brightness tem- 
perature at a given frequency to that predicted from an 
extrapolation of the brightness temperatures measured at 
the three higher frequency points. The disruption caused 
by the reionization signal to the otherwise smooth spec- 
trum is obvious. 
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Fig. 4. The "trend ratio" as described in the text, computed 
for the curves used in Fig. 3, plotted against observing fre- 
quency. 



Clearly, both the magnitude and sharpness of the 
reionization step are critical parameters. As mentioned 
at the beginning of Sect. 3, the magnitude of the step 
will be greater in a low density Universe. The sharpness 
depends on the dumpiness of the Universe at the reioniza- 
tion epoch and the formation rate of the ionizing sources, 
and may well be greater than indicated by curves (b) and 
(c). Simulations of the evolution of structure can give some 
indication of these parameters, but ultimately they have 
to be determined from observations. 

In summary, it appears that the reionization step may 
be detectable in principle, even in the presence of contam- 
ination due to the galactic and extragalactic foregrounds. 
How easily it can be detected, and which method of anal- 
ysis is best, will depend on how well the different con- 
taminating components can be measured, modeled and re- 
moved. For this purpose, new and accurate studies of these 
foreground contaminants will be very important. The re- 
maining limitations are then purely technical, and present 
an interesting challenge. Some of the practicalities are dis- 
cussed below. 

3.3. Observational prospects 
3.3.1. Telescope requirements 

As (1) the signal is present over the entire sky, (2) the 
system temperature is dominated by the galactic and ex- 
tragalactic nonthermal emission, and (3) this is a broad- 
band (rather than line) observation, this experiment can 
in principle be done using already existing radio telescopes 
of moderate aperture. The main requirement is to achieve 
highly accurate relative intensity measurements over a 
considerable frequency range (70-240 MHz). 



Attractive as the use of a small telescope for this mea- 
surement may be, however, there are significant difficul- 
ties. The gain of the telescope, its beamshape and its (dis- 
tant) sidelobes are frequency-dependent, due to match- 
ing, feed illumination, feed support scattering and edge- 
refraction effects. The beamsize increases with decreasing 
frequency, so the regions contributing to the galactic and 
extragalactic emission are different at different frequen- 
cies. This will further complicate the shape of the over- 
all spectrum, although the effect should vary slowly and 
smoothly with frequency, particularly in the directions of 
lowest galactic emission. The distant sidelobes may also 
introduce (varying) signals from regions of stronger galac- 
tic emission. 

The large beam of a small telescope would average over 
the effects of many extragalactic sources, and more sources 
would be included at lower frequencies where the beam is 
larger. Averaging over many sources could conceivably be 
an advantage, as long as a small number of sources did not 
dominate, but with a larger telescope and smaller beam 
it would be easier to observe in regions devoid of strong 
extragalactic sources. Finally, if a strong source such as 
Cas A is to be used as a calibrator (see below), a fairly 
large antenna is required for the signal from this source to 
dominate over the surrounding galactic emission. 

A larger parabolic antenna would solve some of these 
problems, but not all. Some sort of phased array with 
beam-forming circuitry which is still sensitive to emission 
on the largest scales may have several advantages. It would 
be possible to use scaled beams at all frequencies and con- 
trol sidelobes, so that exactly the same area of sky would 
be sampled at all frequencies. New arrays of this type, 
such as the THousand Element Array (THEA) in The 
Netherlands and the Square Kilometer Array (SKA) are 
presently under consideration (van Ardcnnc 1998). 

If terrestrial complications such as radio-frequency in- 
terference were considered insurmountable, one might con- 
sider a small satellite project. Then there would be only 
the galactic and extragalactic foreground contamination 
to contend with, as well as calibration (although solar ra- 
diation may also be a significant problem). The antenna 
size would obviously be limited, but even a deployable 
10m antenna, giving a 10° beam at 180 MHz, may be ad- 
equate. At these frequencies the precision of the antenna 
surface would not be an issue. While space projects are 
not inexpensive, a space radiotelescope such as this may 
find novel uses, as it would be exploring a new domain in 
observational phase space. 

3.3.2. Calibration options 

To be specific we will discuss three different calibration 
strategies, and mention specific problems for each of them 
and how they might be solved. Calibrating the frequency 
response of the telescope/receiver combination is probably 
the crucial issue. This has to be done to better than a 
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few parts in 10 5 over a wide frequency range. Possibilities 
(which could be combined) include internal loads, or the 
use of astronomical sources - a very strong 'featureless' 
continuum source such as Cas A, or the moon. We discuss 
each of these in turn. 

a) Internal loads 

For precise calibration over many frequency channels, 
internal broad-band noise sources of the highest quality 
would undoubtedly be the best option. They would have 
to be extremely stable and very accurately calibrated as a 
function of frequency. The properties of available internal 
calibration sources will have to be examined carefully to 
determine whether this is an option - until now, calibra- 
tion of this accuracy and stability has not been required 
at these frequencies, so this is uncharted territory. Cali- 
bration with an internal load must also take care of any 
frequency dependence of the load-coupling into the signal 
path. 

b) Using Cas A 

The alternative is to use astronomical calibration 
sources. This has the advantage that the astronomical sig- 
nal and the calibration signal both traverse the same path 
through the telescope and receiver system and therefore 
share similar gain-frequency effects. Cas A is the strongest 
source in the sky at 150 MHz, with a flux density of 13,000 
Jy. In order for this source to produce an antenna tem- 
perature well in excess of the local galactic emission (Cas 
A sits on the galactic ridge), a telescope of 25m diame- 
ter or more would be required. The use of such a source 
may introduce a further complication: although its over- 
all power-law spectrum is remarkably straight down to 
a few MHz (Baars et al. 1977), many small components 
with a range of spectral indices are superimposed, mod- 
ifying the detailed overall spectrum to an extent yet to 
be determined (Anderson & Rudnick 1996, and references 
therein). Strong extragalactic sources such as Cygnus A 
have more complex, curved spectra (Baars et al. 1977), 
and so may be less suitable as calibrators. 

Probably the best existing telescope to begin with is 
the GMRT. Cas A would produce an antenna tempera- 
ture of about 4300 K at 150 MHz with the 45m anten- 
nas, totally dominating the sky and receiver noise. The 
GMRT has receivers over a wide range of frequencies, but 
not contiguous over the desired range from 70-240 MHz. 
Total power measurements, accompanied by autocorrela- 
tion spectra of the input signals to identify data affected 
by narrow band interference, are required. Digital spec- 
trometers with at least 20-40 MHz bandwidth and several 
thousand spectral channels would be needed. 

c) Using the moon 

An advantage of the moon over Cas A is that it moves, 
and a true differencing experiment could therefore be done 
at exactly the same position(s) in the sky, and in the same 
alt-az position(s), but of course not at the same time. All 



stable contaminating effects (galaxy, discrete sources, side- 
lobes etc) will to first and perhaps second order be identi- 
cal and will cancel, leaving only time-dependent concerns 
- overall stability and radio-frequency interference. A true 
differencing experiment such as this may be very attrac- 
tive in coping with the contribution to the system tem- 
perature coming in via distant sidelobes. It will be hard 
to limit this contribution to less than 50 K at the low 
frequencies of this experiment. This component may have 
significant frequency structure which could be difficult to 
model. This may be a great advantage of using the moon 
compared to an experiment using Cas A. 

The brightness temperature of the moon is about 220 
K in this frequency range (Keihm & Langseth 1975); it 
varies slowly with lunar phase. This is close to the mini- 
mum sky temperature at 150 MHz, so if the moon fills the 
beam it will not increase the antenna signal by much, if at 
all. The fact that the moon completely blocks the emis- 
sion behind it (galactic and extragalactic) can be very 
useful, and could help solve the problem of the frequency- 
dependent beamsize. The moon is not a very bright cali- 
bration source so it takes a large telescope to avoid dilut- 
ing the signal too much. For the moon to fill the beam, 
we require a telescope of about 200-400 m diameter at 
the frequencies that we are considering. Currently only 
Arecibo would be able to provide a sufficiently narrow 
beam. A specially designed low-frequency phased-array, 
such as those mentioned in Sect. 3.3.1, may be required. 

An interesting aspect of a moon calibration experi- 
ment is that the detection of the signal could in principle 
be done interfcromctrically. Because the moon occults the 
background signal it introduces spatial structure in the 
global edge signal. Depending on the frequency and the in- 
tensity of the galactic background, the 220 K moon could 
be a colder or warmer spot in the sky. But the intensity 
of this negative or positive source contains the frequency- 
dependent global background temperature. 

Of course, the motion of the moon also introduces 
smearing effects, which would probably be advantageous. 
If the 'edge' signal is universal, as expected, we need not 
worry that the moon moves by approximately its diame- 
ter every hour. Both the sought-after signal and the side- 
lobes will be averaged. Since we need to integrate over 
many (tens of) hours by tracking the moon we are effec- 
tively sensitive only to a large-scale signal. The possibil- 
ity that emission from the moon, which includes a time- 
variable contribution from scattered solar radiation, may 
have weak frequency structure will have to be investigated. 

Alternatively (especially if smaller telescopes are used) 
even the sun might conceivably serve for calibration pur- 
poses. Although the solar spectrum has considerable fre- 
quency structure during bursts, the quiet sun is less likely 
to have permanent spectral features at 70-240 MHz, since 
we are looking well into the "quiet" spectrum of the so- 
lar chromosphere superimposed on the Raylcigh- Jeans tail 
of its 6000 K atmosphere. Several issues will have to be 
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further studied before deciding on the optimal calibration 
strategy. 

3.3.3. Terrestrial interference and complications 

Radio-frequency interference (RFI) is of course a major 
issue for observations at the frequencies of interest. The 
World Radio Conference had allocated only one protected 
radio astronomy band in this regime, at 150-153 MHz, 
which itself is already severely affected by interference 
(Spoelstra 1998). This region of the radio spectrum is sub- 
ject to local, regional and world-wide interference. Some 
of these are stationary sources of interference such as TV 
stations, many other in the 100-200 MHz regime are of a 
mobile nature such as aeronautical and military transmit- 
ters and beacons. Fortunately, many of these sources are 
strongly time-variable and confined to relatively narrow 
bands. The 88-108 MHz band is allocated to FM radio 
and may be completely inaccessible for radio observations 
of the required accuracy. This means that if the Hi 'edge' 
signal is at redshift greater than 10 (i.e. at frequencies 
less than 130 MHz) there may not be a reliable spectral 
baseline left. 

It is interesting to consider the role that could be 
played by interferometers in this experiment. Such arrays 
can, with many independent monitors of the total power, 
provide an excellent discriminator (filter) for local or ex- 
ternal RFI. Because RFI generally does not come from the 
direction of the signal its delay and fringe rate will also be 
different. This then provides a way of detecting and pos- 
sibly eliminating RFI, so using telescopes that are part 
of an intcrfcromctric array may yield powerful diagnostic 
capabilities for weak RFI. 

Ionospheric effects are also important at these frequen- 
cies. The ionosphere not only refracts the signal by a vary- 
ing amount (easily 1-2 arcminutes) but there is also refrac- 
tive focussing and scintillation. These scale with frequency 
in a quadratic way. However, the ionosphere changes daily, 
so one can always select the best conditions. The effects 
would be of greatest concern for the observations of the 
discrete calibration sources, and should not be important 
for low-resolution observations of the frequency depen- 
dence of large areas of sky. 

4. Detectability of a Lyman reionization edge 

In this section we discuss the possibility of detecting the 
reionization edge via the Lyman lines. We first discuss the 
nature of this optical/IR signal, next its observability in 
general, and then we conclude by discussing an experiment 
with existing equipment. 

The rate of hydrogen recombinations peaks at z ~ z- lon , 
and this is the source of the background we seek. Because 
we want to observe a spectral signature in the all-sky back- 
ground, we are only concerned with how energy is redis- 
tributed among different spectral lines. Baltz et al. (1998) 



showed that intergalactic absorption brightens some lines 
while it dims others. At z > z- lon , the optical depth in the 
Lyman series is very high. All Lyman series lines except 
Lya are absorbed immediately and redistributed. This 
makes Lya and the Balmer series significantly brighter, 
because they receive the energy of the Lyman series. Once 
reionization occurs, the Lyman series becomes optically 
thin and the brightness of Lya and the Balmer series drops 
considerably, producing sharp features in the background 
spectrum. 

If reionization occurs abruptly, such a sharp edge may 
be observable as a global signal over the entire sky due to 
Lya emission at z~5-20 (A ~0.7-2.6/im). This can be seen 
from the models of Baltz et al. (1998) shown in Fig. 5. As 
we look redward, crossing Lya around ~l/xm during the 
reionization event, the brightness suddenly increases due 
to the absorption and redistribution of the higher Lyman 
series lines. Haiman et al. (1997) calculated a similar effect 
in which Lya is brightened by radiative excitations, but it 
appears that recombination gives a significantly stronger 
signal. Baltz et al. (1998) and Longair (1995) estimate 
the strength of the edge to be in the range J„(Ly) ~ 0.3 — 
3 x 10~ 23 erg/cm 2 /s/Hz/sr (these estimates are probably 
uncertain by £ 0.5 dex). 

Further recombination features are also present in the 
spectra in Fig. 5. The broad humps (most noticeably that 
blueward of Lya) are due to blends of lines from the Ly- 
man and Balmer series (Baltz et al. 1998; see also Stancil 
et al. 1998 for other possible transitions with lower prob- 
ability). The positions of these humps are of course fixed 
relative to the sharp features. Their shapes are determined 
by the evolution of the reionization/recombination history. 

4-1. Fundamental sensitivity limits 

The sky brightness at these wavelengths is dominated 
by the zodiacal light. The Diffuse Galactic Light (DGL), 
and the Extragalactic Background Light (EBL) also con- 
tribute, but at a lower level. These "foreground" emissions 
determine the fundamental sensitivity limits for this ob- 
servation. 

The dominant foreground signal is the zodiacal "back- 
ground", which at high galactic and ecliptic latitudes is 
well determined by HST in the WFPC2 B, V, I bands 
at -23.9, 23.00, and 22.45 arcsec -2 , respectively (Wind- 
horst et al. 1992, 1994b, 1998), and £20.9 and £19.6 
mag arcsec -2 in the Near-Infrared Camera/Multi Object 
Spectrograph (NICMOS) ncar-IR J and H bands, respec- 
tively (Thompson et al. 1999; Windhorst & Wadding- 
ton 1999). The latter are upper limits to the IR zodi- 
acal sky, since there is some remaining uncertainty in 
the global NICMOS dark-current. At 0.55yum, the zodi- 
acal sky as measured from low earth orbit corresponds to 
J„ — 6.8 x 10~ 19 erg/cm 2 /s/Hz/sr, and at 1/xm it corre- 
sponds to 3.1x 10~ 18 erg/cm 2 /s/Hz/sr. 
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Fig. 5. Predicted spectra of the Lyman reionization edge 
from Baltz et al. (1998). Models A-D were generated from 
CDM simulations including radiative transfer, star-formation, 
feedback, evolution and chemistry of the gas, and assume 
&m = 0.35, = 0.65, and h — 0.7. The sequence A-D has 
decreasing values of Qb (0.055-0.03), simulation box-size (3-2 
h -1 Mpc), mass resolution (10 7 '°-10 6 '° M Q ), and number of 
particles. The amplitude prediction in the models is probably 
uncertain by <; 0.5 dex. The most noticeable hydrogen edges, 
Lya and Ha, are indicated (for model D only). The hump blue- 
ward of Lya is caused by higher atomic transitions. The dotted 
curves indicate the in-principle 3-<r limits expected from a day 
of integration with HST/STIS, HST/WFC3, NGST, and a long 
integration on a similar spacecraft at 3 AU from the sun. The 
uncertainty in these limits are about ~0.3 dex. 



At high galactic latitudes and in regions of low Nh, 
the DGL contributes less than 26.0 mag arcsec -2 at 
UV-optical wavelengths, or ^ 6 x 10 -20 erg/cm 2 /s/Hz/sr 
(Henry 1991, O'Connell et al. 1992). 

The EBL surface brightness is estimated to be 
IJJ24-25 mag/arcsec 2 , since the galaxy counts are now 
known to permanently converge with a sub-critical mag- 
nitude slope a ^0.4 for BJ>25 mag (Metcalfe et al. 
1995) and for L>24 mag (Williams et al. 1996). The 
best estimate of the extragalactic background is about 
1.7x 10 -20 erg/cm 2 /s/Hz/sr or ^27.3 mag arcsec -2 at 
0.55/im, and 6.7x 10 -20 erg/cm 2 /s/Hz/sr or ^25.1 mag 
arcsec -2 at lfim (Dwek et al. 1998; Hauser et al. 1998) 

It is assumed here that the high-redshift intergalactic 
medium is not obscured by dust and neutral clouds at 
redshifts somewhat lower than Zi on . Comparisons of radio 
and optical samples of high-redshift quasars (Shaver et al. 
1996, 1999) and the discovery of Lya emitters at z~5-6 
(Dcy et al. 1998, Weymann et al. 1998) support this as- 



sumption. In any case, while foreground dust could dimin- 
ish the observed amplitude of any reionization signal, it 
would not significantly affect its global spectral signature 
(c.f., Scaton 1979). 

Taking the above foreground emissions into account, 
the horizontal lines in Fig. 5 show the sensitivity expected 
from a one-day integration with HST/STIS (Space Tele- 
scope Imaging Spectrograph), HST/WFC3 (Wide Field 
Camera 3), NGST (the Next Generation Space Telescope), 
and a long integration on a similar spacecraft at 3 AU from 
the sun where the zodiacal emission would be reduced by 
two orders of magnitude. The reionization signal is within 
reach, at least for the deep space mission. 

We note that in ground-based direct CCD-imaging lo- 
cal sky-removal to a precision of <~ 10 -3 is obtained rou- 
tinely, and that it has been obtained to a precision bet- 
ter than ~ 10 -4 when using very careful calibrations of 
the kind we propose below (e.g. Tyson 1988). There is no 
fundamental reason why a CCD spectrograph outside the 
Earth's atmosphere should not be able to achieve these 
limits in spectroscopic mode, although this has yet to be 
demonstrated. 

The sensitivity requirement may be eased by taking 
advantage of the considerable spectral structure present 
in the curves shown in Fig. 5. Rather than looking only 
for the Lya and Ha steps, it may be possible to search 
for the overall pattern in these curves by cross-correlating 
the observed spectrum with a theoretical template (e.g. 
Baltz et al. 1998), as is often done in determining redshifts 
of elliptical galaxies. This could ease the sensitivity re- 
quirement by a considerable factor (which depends on the 
actual structure in the reionization signal). In this case, 
however, the observable would only be a cross-correlation 
peak rather than the actual spectral features, and such a 
technique would require excellent flat-fielding and calibra- 
tion over a wide wavelength range. 

4-2. Contamination by galactic and extragalactic 
foregrounds 

Any spectral features or variations in the foreground con- 
taminants could conceivably mask the reionization signal, 
but it is likely that such contamination can be adequately 
dealt with. Deep ground-based imaging of each field ob- 
served would identify any large-scale low-surface bright- 
ness objects or structures to ^ 29-30 mag arcsec -2 (e.g. 
Tyson 1988) that can be excluded afterwards. The ob- 
served sky background spectrum can then be corrected 
for the diffuse galactic and extragalactic foreground emis- 
sions as follows: 

(a) The zodiacal (~solar) spectrum is the main con- 
tributor to the variations in the spectral baseline, and is 
mostly featureless over the wavelength range of interest, 
except for perhaps the solar Ca feature at 0.89yum. The so- 
lar spectrum is known to very high accuracy. The zodiacal 
dust particles with sizes <l^tm do not significantly change 
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the high frequency part of the zodiacal spectrum around 
1/zm, but they slightly modify the global solar spectral 
gradient due to non-grey scattering. Such changes are not 
expected to be sudden with wavelength. A slightly red- 
dened solar spectrum can be subtracted from the final 
observed sky spectrum after appropriate scaling to fit the 
HST broad-band zodiacal sky-measurements. 

(b) Possible high spatial frequency structures in the 
galactic DGL at very low surface brightness levels may be 
of some concern, but only at a level of >27 mag/arcsec 2 , 
as has been seen in deep ground-based imaging at 
low galactic latitudes (Szomoru & Guhathakurta 1998; 
Guhathakurta & Tyson 1989). At high galactic latitudes, 
if the observed fields can be selected to have Eb-v & 0.01 
and/or low Njj, the DGL and variations therein should be 
very small, as shown by deep ground-based images (e.g., 
Tyson 1988) and deep HST images such as the Hubble 
Deep Field (Williams et al. 1996). Fields near potentially 
contaminating stars or galactic cirrus would have to be 
avoided. Since the DGL likely bears the spectral signa- 
tures of early-type stars in our Galaxy, we do not expect 
a large spectral feature around 1/xm from the combined 
stellar population that causes the reflected DGL. In any 
case, such a signal can be subtracted from the residual sky 
spectrum at the B~27 mag arcsec -2 level by subtracting 
a scaled OB star template, as for the zodiacal spectrum. 

(c) Finally, the integrated background from faint fore- 
ground galaxies (EBL) should be featureless, given their 
large redshift range (cf. Driver et al. 1998). Therefore no 
global template for the EBL from faint galaxies at z <, 5 
should have to be subtracted from the collapsed spectrum 
to find a Ly-edge. 

4-3. Observational prospects 

4.3.1. Observing the Lyman edge with space-based 
instruments 

A near-infrared instrument in space would be essential for 
this experiment, because of the low sky background (in 
space the zodiacal sky in the H-band is 6-7 mag darker 
than from the ground (Thompson et al. 1999; Windhorst 
& Waddington 1999) and complete absence of terrestrial 
OH lines. We note that a recent deep long-slit search with 
Keck provided quite useful upper limits to the ionizing 
background J„o at levels 2 x 10~ 21 erg/cm 2 /s/Hz/sr, but 
at lower redshifts (2.7^z^4.1; Bunker et al. 1998). In our 
case we are concerned with redshifts z^5, and so we must 
do the experiment from space to avoid the strong OH- 
forest. 

Ideally, one would do this experiment on a faster tele- 
scope with a larger aperture and lower sky than HST, such 
as the NGST, to be launched sometime after 2007. Im- 
proved infrared cameras with spectroscopic multiplexing 
capability, such as those under consideration for WFC3 
and being designed for NGST, are also essential. They will 



have larger IR detectors with much better noise character- 
istics than HST/NICMOS, making possible a substantial 
gain in sensitivity for this experiment. The NGST will op- 
erate over the range A =0.5-10/im, and will have a large 
spectroscopic multiplexing capability provided by an inte- 
gral field or multi-object spectrograph with some 40,000 
spectral pixels, One would search for the isotropic step- like 
spectral features in the collapsed 1-dimcnsional spectrum 
of the sky. 

Alternatively, an infrared space telescope dedicated to 
studying the CIBR should be able to detect the Lyman 
edge signals. Two versions, EGBIRT and DESIRE, have 
already been proposed (Mather & Bcichman 1996). These 
would observe at 3 AU from the sun, reducing the zodiacal 
foreground by two orders of magnitude. 

The sensitivities attainable with such space-based fa- 
cilities are shown in Fig. 5. As indicated above, the level of 
the reionization signal can be reached, but it may require 
the "deep-space" mission with its much lower background 
and long dedicated integrations. In the meantime, an ex- 
ploratory experiment can be done with HST, as described 
below. 

4.3.2. An experiment with the HST 

In view of the large uncertainties, both in the predicted 
amplitude of the reionization signal and in the attainable 
sensitivities, it is worth considering preliminary experi- 
ments using the currently existing Space Telescope Imag- 
ing Spectrograph (STIS). Such experiments would also be 
exploring new observational domains, and unexpected dis- 
coveries may result. For example, in addition to the hy- 
drogen recombination signal shown in Fig, 5, there may 
also be a broader and possibly stronger spectral feature 
due to the continuum emission from the ionizing sources 
near Zi on as shown in fig. 11 of Gnedin & Ostriker (1997). 
Such an experiment would in any case provide important 
experience in making observations of this kind. 

The STIS has a wide variety of long slits, permitting 
a spectral image of the sky to be obtained. To obtain 
sufficient surface brightness sensitivity, a long observa- 
tion with the widest available STIS long slit of 52" x 2.0" 
should be made. One would use the longest-wavelength 
low-resolution STIS grating G750L to make a deep spec- 
trum of the sky. The observation would be limited by the 
STIS CCD to <> 1.03/xm, or z <, 7.5 for Ly-a. 

This HST observation could be done in parallel mode, 
at no extra cost in primary HST observing time. As 
the primary exposures would most likely be taken with 
WFPC2 and therefore be dithered, the STIS exposures 
would also be dithered along the slit, resulting in better 
spectral flatfields. Several independent fields could be ob- 
served in order to identify and remove "spurious" features 
from foreground objects. One or two fields should be taken 
at lower ecliptic latitudes - or higher zodiacal foreground 
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- so that the zodiacal spectrum could be removed better 
by comparison with the higher latitude observations. 

A project such as this relies on the long temporal sta- 
bility and excellent flat-field properties of ST1S, and would 
push the ST1S calibration to its limits. Contemporaneous 
calibrations would be essential to reduce the systematic 
errors in the spectra to an absolute minimum. The entire 
bright side of each orbit should be used to make the max- 
imum number of calibrations between the parallels taken 
during the dark side of each orbit, without interfering with 
the primary observations: bias frames, dark frames, earth 
flats (done between the two darks in the middle of the 
bright side of the orbit, so that HST would be looking 
down for earth flats), and internal tungsten flats. 

The experiment would be read-noise limited, since one 
cannot expose for long periods and still reject cosmic rays. 
Sensitivity could be gained by on-chip binning, as the ob- 
servations are not dark-current limited (Baum et al. 1998). 
Fringe removal could be done with contemporaneous spec- 
tral flats (Baum et al. 1998). Following cosmic-ray removal 
(c.f., Windhorst et al. 1994a), the best possible bias, darks, 
and flats would have to be produced. The end goal would 
be a spectral image that is globally flat-fielded to within 
< lCT 3 x sky. 

The following sensitivities can be reached in princi- 
ple with the STIS/CCD/G750L combination, according 
to the STIS Manual (Baum et al. 1998) and the STIS 
exposure time calculator. With 36 orbits (total time 94 
ksec) a la surface brightness sensitivity of 1.4 x 10 -20 
erg/cm 2 /s/Hz/sr can be achieved per pixel, or a 3<7 sen- 
sitivity of 1.3 x 10~ 21 erg/cm 2 /s/Hz/sr if the sky is flat to 
well within ^ 10 -3 , so that one could collapse all 1000 
spatial pixels along the slit into a one-dimensional STIS 
sky spectrum. The approximate 3-cr limit expected from 
a 36 orbit HST/STIS integration with the grating G750L 
and the 2" long-slit is shown in Fig. 5. The uncertainty in 
this STIS simulation is about 0.3 dex. 

The question then is whether any significant spectral 
structure survives in the collapsed STIS spectrum after 
subtraction of the best-fit scaled templates of the zodiacal 
and DGL spectra, and after checking for possible coin- 
cidences with any of the weak known Ca features in the 
solar spectrum (or perhaps in the DGL) around 0.89/im. If 
structure is not obvious, large wavelength swaths could be 
summed and moved in a sliding manner to find the most 
significant features. Alternatively the spectrum could be 
cross-correlated with a theoretical template such as the 
curves in Fig. 5, or autocorrelated to find any possible 
steps buried in the data. Such an experiment, pushing 
STIS to its limits and exploring a new observational do- 
main, will be an important first step both technically and 
scientifically. 



5. Conclusions 

The reionization of the Universe may be detectable glob- 
ally in the form of a step in the broadband spectrum of 
the sky. This would be a unique measurement of a basic 
"phase transition" of the Universe. As such it would differ 
fundamentally from previously proposed searches based 
on fluctuations due to individual H I concentrations or ab- 
sorption against continuum sources which may exist at 
still higher redshifts. The step is expected at ~ 70-240 
MHz due to redshiftcd Hi 21 cm line emission and at 
A ~ 0.7— 2.6/im due to hydrogen recombination emission. 
The expected amplitude is well above fundamental de- 
tection limits. The difficulties lie in contamination from 
galactic and extragalactic foregrounds, calibration, and 
terrestrial interference. 

It appears that complexities due to the frequency de- 
pendence of the galactic and extragalactic foregrounds 
may not be insurmountable obstacles to identifying and 
measuring the reionization step. If that is the case, the 
experiment is possible in principle, and the remaining hur- 
dles are purely technical. The challenge may be justified 
not only by the potential importance of this measurement, 
but also by the possibility of opening up new domains in 
observational parameter space. 
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